Interface and contact line motion in a two phase fluid under shear flow 
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We use a coarse grained description to study the steady state interfacial configuration of a two phase 
fluid under steady shear. Dissipative relaxation of the order parameter leads to interfacial slip at 
the contact line, even with no-slip bound ary con ditions on the fluid velocity. This relaxation occurs 
within a characteristic length scale lo = y/^D/Vo, with £ the (microscopic) interfacial thickness, D an 
order parameter diffusivity, and Vo the boundary velocity. The steady state interfacial configuration 
is shown to satisfy a scaling form involving the ratio lo/L, where L is the width of the fluid layer, 
for a passive interface, and the capillary number as well for an active interface. 



We reexamine a classical hydrodynamic problem in 
which a two-phase fluid is placed in a uniform shear 
flow by displacing two infinite, parallel planar boundaries 
at constant relative velocity. Instead of considering the 
case of two immiscible fluids separated by an infinitely 
sharp interface, our analysis involves two coexisting fluid 
phases. The interface, which is diffuse even without 
shear, is stretched by the flow, but not indefinitely as 
interfacial slip relieves the increasing shear stress at the 
three phase junction. A mesoscopic description of the 
two-phase fluid is introduced to show that diffusive re- 
laxation of the order parameter (e.g., concentration) nec- 
essarily contributes to macroscopic interfacial slip. Since 
in typical situations mutual diffusion occurs on lengths 
large compared to microscopic (molecular) lengths, it is 
not necessary to introduce microscopic slip of the fluid 
velocity near the contact lines. Finally, we also show that 
in the steady state, the interfacial configuration satisfies 
a scaling form, provided that the interfacial curvature is 
small compared to its inverse thickness. 

The classical analysis of the contact line profile and 
motion in partially wetting fluids at low capillary num- 
bers is due to Cox and Dussan and Davis g]. A 
matched asymptotic expansion describes both the inter- 
facial configuration and flow field by considering three 
separate regions. Close to the boundary, the fluid is al- 
lowed to slip by invoking a phenomenological relation 
between the slip velocity and the local shear stress. This 
relation can also be viewed as introducing a slip length 
of microscopic dimensions and an effective slip boundary 
condition for the velocity. Far from the boundaries, capil- 
lary forces are balanced by viscous stresses. An interme- 
diate region allows matching between the two. Reviews 
of this treatment have been given in . Early experi- 
ments by Dussan J5| confirmed this picture. More recent 
experiments by Lichter || have addressed the range of 
capillary numbers within which this solution is a good 
approximation near the contact line. 

If the fluid in question is a binary at coexistence, 



interfacial deformation induces concentration gradients 
through a Gibbs-Thomson relation Q| . The resulting dif- 
fusive transport leads to effective interfacial slip at the 
contact line even when the velocity strictly satisfies a no- 
slip boundary condition. Increasing interfacial curvature 
leads to increasing diffusive transport which eventually 
balances advection by the flow. A steady state can be 
maintained. Within a coarse-grained or mesoscopic de- 
scription of the two phase fluid H, this balance occurs 
within a characteristic length scale lo, which is the ge- 
ometric mean of the interfacial thickness and the diffu- 
sion length, and hence is a small quantity although not of 
molecular dimensions. Dissipation within this scale effec- 
tively prevents the formation of stress singularities at the 
contact line. A mesoscopic description of contact line mo- 
tion has already been given || , by using the same model 
equations which we describe below. These authors, how- 
ever, focused on the limit /o — ^ (in our notation), which 
corresponds to effective immiscibility on all scales. The 
effect that we discuss, however, requires consideration of 
distances to the boundary small compared to ^o- 

We consider an incompressible binary mixture in two 
spatial dimensions confined between two planar bound- 
aries of infinite extent located at x = ±L/2. The 
bounding planes are displaced with constant velocity 
±Vb/2 y. We study the steady state interface config- 
uration in terms of Voj L, and the physical parameters of 
the fluid. Within a coarse grained description of a binary 
fluid 1 10 13, the order parameter <f>, e.g., concentration, 
satisfies that <p = <p + > for one phase, and (j> = <f>- < 
for the other. The nominal position of the interface is 
defined as the locus (f> — 0. The temporal evolution of 
the order parameter is governed by 



dt 



• v • = AfV> , 



(1) 



where v is the fluid velocity, the chemical potential 
= SF/Scf), with F\<t>] the coarse grained free energy 
functional, and 5/5<f> stands for variational diffcrentia- 
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tion. M is a mobility coefficient assumed to be constant. 
The velocity v obeys a modified Navier-Stokes equation, 



'dv „ 
,M W +v.Vv 



-Vp + /xV0 + 7]V 2 V 



(2) 



where p and 77 are the density and shear viscosity of 
the fluid, assumed to be independent of </>, and p is the 
pressure. The second term on the right hand side of 
Eq. (0) incorporates the effects of capillarity. We choose 
the standard free energy F[4>] = J dV (f \V(/)\ 2 + /(</>)), 
with K constant. The first term represents the excess 
free energy due to spatial inhomogeneities of <j>, and f((f>) 
is the local part of the free energy density. Generaliza- 
tions involving, say, higher gradients or forcing due to 
short-ranged wall interactions can be incorporated. 

Several macroscopic physical quantities follow from 
this coarse grained description JijJt]]. The surface ten- 
sion is given by a — K [defio/dy] dy, where (po is 
the equilibrium order parameter profile at coexistence 
(we take fi(<f>o) = in what follows). Furthermore, the 
chemical potential at a gently curved interface with ra- 
dius of curvature R c , relative to its value at coexistence, 
is (i = —cr/R c A<f), where A<j) = <f> + — <f>_ is the miscibility 
gap. Another important quantity is the (finite) interfa- 
cial thickness which is proportional to £ = v*Kx where 
X -1 = (<9 2 //90 2 )^ ± . Since we assume that the fluid 
is not too close to criticality, £ is of molecular dimen- 
sions. Finally, we focus exclusively on moderate shear 
rates, such that £/i? c <C 1, a limit commonly realized in 
practice. 

We begin by illustrating the effect of the shear on the 
interface through dimensional analysis. If the interfa- 
cial mean curvature is l/l, then the local deviation of 
fi from coexistence is proportional to l/l. We therefore 
choose the following dimensionless (barred) quantities: 
fi = p, <j) = (A(f>)4>, and v = V v. In steady 

state, Pew ■ V<f> = V 2 fl, where Pe = V A(f) 2 l 2 /Ma. On 
small scales (I ~ £), Pe becomes very small and hence 
V 2 /2 = defines the equilibrium order parameter pro- 
file in the absence of shear. Diffusion and miscibility 
are important. If, on the other hand, I ~ L — > 00, then 
V^>-v = 0, corresponding either to a homogeneous config- 
uration, or to an interface parallel to the flow. Miscibility 
effects are not important at this scale. There exists, how- 
ever, an intermediate length scale Iq = y 'Ma /VqAc/) 2 for 
which Pe = 1. Advection now balances order parameter 
diffusion along the interface driven by chemical poten- 
tial gradients; this is the length scale of interest. Given 
the definition of the interfacial thickness £, we can write 
lo oc W^DJVq , where D = M/x is a diffusion coefficient. 
Hence this new length scale is the geometric mean of the 
interfacial thickness and a diffusion length D/Vq. For a 
typical system, D ~ 10~ 5 cm 2 / s, £ ~ 10 A, and for a wall 
velocity Vq ~ 1 cm/s, Iq w 10 nm. 

We next present an approximate analytic treatment 
that is valid for length scales larger than £. We find 
a scaling solution for the steady state interfacial profile 



for both passive interfaces (neglecting capillarity induced 
backflow, i.e., model B in the critical dynamics lexicon) 
and active interfaces. We then present a numerical so- 
lution of the full set of governing equations to validate 
the scaling forms derived, at least within the range of 
parameters which is numerically accessible. 

For a passive interface, the velocity field is given by 
v s (r) = xVo/L y. Equation (|l|) in steady state is written 
in integral form, 



M/i(r) = J dr'G(r, r') [v • V<; 



(3) 



where G(r,r') is the Green's function of the Laplacian 
operator satisfying Neumann boundary conditions at x = 
±L/2 |15| . In two dimensions the Green's function is 
dimensionless. Let the position of the nominal interface 
be 17 = yi(x)y + xk, where x and y are unit vectors. 
When radii of curvature are large compared to £, one 
can make the standard approximation, ~ (A(f>)5(r — 
r/)n, and = 6±. jl6| Using now standard methods 
(see, e.g., Ref. ]7|,|l6||), one multiplies Eq. (||) by d<f>/dg 
(g = r • V^/|V0|) and integrates g across the interface 
to derive an equation that depends on the coordinates of 
the interface alone. The result is 



R(X) 



1/2 



1/2 



dX'X'G(Hi, Hi) , 



(4) 



where we have now scaled lengths by L (R = R c /L, 
X = x/L, Y = y/L), and Ri(AT) = rj(x)/L is the scaled 
location of the interface. This equation is an integro- 
differential equation for the interface coordinate Yj(X). 
Boundary conditions need to be imposed at a distance of 
order £ from the boundary. At this scale the shear is neg- 
ligible so that the equilibrium order parameter profile is 
adequate. For unforced conditions, this leads to Y! = 0, 
although other situations can be accommodated. (l7| Ac- 
cording to Eq. (|^) , the curvature of the interface satisfies 
the following scaling relation 



R- 



l (X) = h((lo/L) 2 , 



X 



(5) 



Hence, the interface configuration is solely determined 
by the dimensionless number (Iq/L) 2 , and by the equi- 
librium boundary conditions on the order parameter at 
x = ±L/2. 

Since Eq. (^) is expected to be valid only on length 
scales large compared to £, we have numerically solved 
Eq. ([[]) to verify the above scaling form, as well as to de- 
termine the approach to this scaling form as io/C 1- 
The details of the numerical algorithm have been de- 
scribed elsewhere [jlCj. Briefly, the fluid is in a rectan- 
gular domain -L/2 < x < L/2 and -3L/2 < y < 3L/2, 
with boundary conditions n • V</> = and n • V(V 2 0) = 
at x = ±L/2, and <f> = <j>± at y = ±3L/2. (The corre- 
sponding equilibrium contact angle is 7r/2; the boundary 
conditions also ensure no flux of <f> through the lateral 
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walls.) The interface is initially located at y = 0, and 
we integrate Eq. (|l|) forward in time until a stationary 
configuration is reached. The local part of free energy 
density is the standard double well potential []l(j, and 
a square grid has been used of side no larger than £/3 
so that the relative error in the location of the interface 
is 5% or less. Fig. [j] shows the interfacial configuration 
for a range of Iq/^. The curves tend to superpose as /o/£ 
becomes larger for fixed lo/ L, supporting the scaling pre- 
diction. 

For the range of parameters considered, the interfacial 
displacement away from planarity is small. Therefore 
we also obtain a perturbative solution of Eq. (Q). Direct 
substitution of the Green's function |l5| into Eq. (ji|) leads 
to, 



Yi(X) 



L 



2 oo 



at lowest order in Yj. Note that the term corresponding 
to n = 1 dominates so that the correct dimensionless 
parameter for the expansion should be (L/lo) 2 ir~ 5 . We 
next define S = —1/2 + X and find the curvature of the 
interface near the boundary, R~ 1 (—l/2 + 6) = —R^, 1 + 
±(L/l fx [(3-2m(7r/2))6 2 -2<5 2 ln<5] +... , where 
terms 0((L/lo) 4 ) and 0(5 A ) have been neglected. The 
constant R w « (0A26/2tt)(L/I o ) 2 ■ Note that there is a 
weak singularity as 6 — > (a divergent second derivative) , 
that can be traced back to the factor (2n— 1)~ 5 in Eq.(f3|). 
The perturbative solution for Yi is also shown in Fig. 1. 
The singularity in the second derivative of the curvature 
is much too weak to be observed directly. 

We next consider an active interface and allow for the 
additional flow induced by order parameter variations. 
Boundary conditions for Eq. (||) are no-slip at the solid 
boundaries, v = ±Vb/2y at x = ±L/2, and a far field 
velocity approaching v s . The steady state configuration 
of fi is again given by Eq. (||), but v needs to be de- 
termined self-consistently. Under the same assumptions 
used to derive the interface equation in the case of a pas- 
sive interface and neglecting the inertial term, one finds 



6, 



1 

RjX) ~ 

1 dS 
Ca'dX' 



lo 
L 

1/2 
-1/2 



1/2 



dX'G(R/,R' 7 ) [-X'+ 

-1/2 



(7) 



Repeated indices imply summation, and TJy is the Os- 
een tensor JT8| . A capillary number has been defined as 
Ca = rjVo/a, — m(X') is the ith component of the 
unit normal to the interface at X', and dS{X) / dX = 
a/1 + (dYj(X)/dX) 2 . This interface equation, in turn, 
suggests the following scaling form, 



R 



- 1 (X) = h(Ca,{h/Lf 



.X 



(8) 



In the limit Ca — > oo, the passive case is recovered 
(Eq. ®). 

We have also verified this scaling form by direct numer- 
ical solution of Eqs. (|l|) and (|J) subject to the boundary 
conditions given: no-slip at x = ±L/2, and v = v s at 
y = ±3L/2. Figure || shows the results obtained for 
fixed Iq/L and Ca and a range of values of 2o/£- Again, 
these results are consistent with the scaling forms de- 
rived. Note that the values of Iq/L, and lo/£ are the 
same in both Figs. |l| and ||, but the interface deflection 
in Fig. U is much smaller than that in Fig. This is a 
consequence of the additional dissipation mechanism in 
the fluid and the induced flow. The small values of Y} in 
the numerical solution suggest that Yj is proportional to 
Vb, as can be seen from the inset in Fig. 2. 

We now compare our results with previous studies 
of immiscible fluids based on a macroscopic descrip- 
tion [0-^1. Equation (||) can be written in the follow- 
ing scaling form l /R c (x/lo) — f2(Ca, Iq/L, x/l ). Now 
choose the origin such that the walls are located at x = 0, 
and L. In the limit x/L « 1 but x/Iq 3> 1, one expects 
that the dependence of the curvature on both L and Iq 
is weak. If the leading behavior is such that the curva- 
ture is independent of both L and ^o, then the scaling 
function reduces to /2 ~ fs(Ca)lo/x. Using the fact that 
fz{Ca) ~ A Ca jl9|], this relation can be integrated to 
yield ^p- = A Q Ca In ^- + B where A and B are con- 
stant, and xo/lo 1. This "outer solution", valid for 
small Ca and x/Iq ^> 1, has the same functional form as 
the immiscible case [0. There, the expansion parameter 
in the matched asymptotic expansion is e = s/L, where s 
is constant and of molecular dimensions. In the present 
case, however, the inner region is determined by lo, which 
is a function of the boundary velocity. Therefore a similar 
matched asymptotic expansion to derive the macroscopic 
interfacial profile for a two phase fluid would lead to a 
qualitatively different behavior. In particular, note that 
as Vb becomes small, nonlocal diffusive coupling between 
distant regions of the interface becomes important. We 
note finally that our conclusions are based on the scaling 
form derived from Eq. (^), which is valid only for x £, 
although it has also been numerically validated for the 
range of parameters that we can explore computationally. 

To summarize, we have studied the steady state inter- 
facial configuration in a partially miscible two-phase fluid 
under steady shear. Including mutual diffusion and mis- 
cibility naturally introduces a cutoff, physical in nature, 
which eliminates stress singularities. An approximate 
integro-differential equation for the steady state interface 
configuration has been derived, valid at distances larger 
than the interfacial thickness £. For a passive interface, 
the configuration satisfies a scaling form that depends 
only on the ratio Iq/L, where Iq = y/^D/Vo. We then 
argue that the length scale lo, which is much larger than 
molecular dimensions but not necessarily macroscopic, 
defines the scale over which effective interface slip occurs 
purely due to diffusive relaxation of the order parameter. 
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Similar results are obtained for an active interface, with 
the scaling function now depending on both Iq/L and 
the capillary number Ca = r/Vo/cr. Numerical solutions 
of the full dynamical equations, Eqs. ([!])- (0), support 
both scaling forms. Our results indicate that the slip 
mechanism of an interface separating coexisting fluids is 
qualitatively different than in the immiscible case. Ex- 
periments involving very small Ca on coexisting fluids 
would be useful to verify our results. 
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To lowest order in Ca, this leads to = — CaX + 

/-i/2 dx ' T vy( x !< x ' hence the curvature is 

proportional to Ca when Ca is small. 
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FIG. 2. Scaled steady state configuration of an active inter- 
face for the values of Zo/C listed, L = 4.04Zo with Ca = 0.053. 
The inset reveals the maximum interface deflection A for dif- 
ferent wall velocities (i.e., different Ca) with fixed L/£ and 
l\Ca. 
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